arXiv: 1509.0022lv2 [cond-mat.stat-mech] 20 Sep 2015 


Classical many-body systems with retarded interactions: 
dynamical irreversibility 


A.Yu. Zakharov^, M.A. Zakharov^ 

Yaroslav-the-Wise Novgorod State University, Veliky Novgorod, 173003, Russia 


Abstract 

The exact equations of motion for microscopic density of classical many-body system with 
account of inter-particle retarded interactions are derived. It is shown that interactions retar¬ 
dation leads to irreversible behaviour of many-body systems. 
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1. introduction 

There is a well known contradiction between reversible equations of classical mechanics 
and irreversible behaviour of many-body systems. In particular, it remains an open question 
about the possibility to describe transition of many-particle system to equilibrium without the 
use of the probabilistic hypotheses. Numerous paradoxes, including the Loschmidt and the 
Zermelo paradoxes [i|, 0 , really prove the incompatibility of deterministic classical mechanics 
and probabilistic approach to mechanics. In order to use the probabilistic concepts in classical 
mechanics, a physical source of random behaviour should be indicated. There is no such 
sources in deterministic classical mechanics. In particular, there is no random sources for 
isolated from external world Gibbsian microcanonical ensemble in statistical mechanics. 

However, it is very doubtful that non-equilibrium many-body system isolated from the 
outside world does not go into an equilibrium state as long as does not interfere a random 
external source. Therefore, we must perform an analysis of unstudied factors that can lead to 
irreversible behavior of many-body systems. There are at least two such factors. 

1. Interactions between (electrically neutral atoms are due to mutual polarization of atoms 
and depend on their polarizabilities [1, 0) S 0 ■ This means that “instant” interaction 
between atoms depend on relatively fast electronic degrees of freedom of atoms. In 
addition, there are around vacuum fluctuations of the electromagnetic field, manifested 
in such well-known effects, as the Casimir effect and the Lamb shift. Therefore, there 
is an additional time-dependent contribution into the interatomic potentials realized 
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through the vacuum fluctuations of the electromagnetic field Q. Thus, in addition to the 
static part of the interatomic potentials there exists fluctuating in time interaction. Note 
that even insignificant in quantity fluctuation correction to the interatomic potentials 
can lead to a qualitative change in the behavior of classical systems due to the instability 
of the dynamical systems trajectories Q- Probably, this effect can be taken into account 
by statistical methods only. 

2. There is the retardation effect of interatomic interactions. Accounting for this effect 
involves the use of systems of functional differential equations, i.e. differential equations 
involving unknown functions for different argument values. The differences between the 
argument values of unknown functions are called argument deviations. The mathemati¬ 
cal formulation of problems for functional differential equations is radically different from 
the problems for ordinary differential equations. Unfortunately, currently the theory of 
functional differential equations is not well developed. More or less effective methods of 
the problem functional differential equations solutions are realized mainly in the control 
theory Q . The hypothesis that the delay interactions can be the cause of irreversibility 
expressed by various authors (see for example). 

The purpose of this paper is to consider the effect of interatomic interactions retardation 
on the dynamics of a classical many-body system. 

Unfortunately, the direct account of retarded interactions in conventional approaches of 
classical mechanics is very troublesome: it is necessary to write a system of ordinary differ¬ 
ential equations of motion for all particles of a many-body system in view of the interactions 
retardation. It is absolutely immense problem. The fact that the usual forms equations of 
classical dynamics describe the motion of each particle of the system. For a many-body system 
such details are redundant. We are interested in the evolution of the particle distribution in 
the space as a whole. 

For this task, we need to describe a many-body system evolution in terms of the local 
microscopic density. Such a form of the equations of motion in the classical many-body 
problem for interactions without delay has been proposed in the papers [3, El 

N r d\< 

= X! = / UTTs (1) 

S = 1 

where Rs(t) is position of s-th particle at time t, N is the number of particles in the system, 
fi(k, t) are the collective coordinates [13 

h(k,t) = (2) 

S 

In this formulation, we do not assume any probabilistic interpretation of n(r, t). Since the 
dynamical equations of classical mechanics determine the evolution of the system uniquely, 
the function n[r,t) is a completely deterministic function. 

The aim of the present paper is derivation of exact equation for many particle classical 
system evolution in terms of the local microscopic density with account of the interatomic 
potentials retardations. This derivation is based on the Newtonian laws only and does not 
based on any probabilistic or statistical assumptions and hypotheses. It is shown that the 
retardation interactions leads to irreversible dynamics of many-particle system. 
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2. Equations of motion for microscopic density of a many-body system 

To obtain the equations of motion for the function n(r, t) we need differentiate this function 
with respect to time and use the Newtonian laws. Since the Newtonian second law contains the 
second derivatives with respect to time, we should also be restricted to the second derivatives 
of n(r, t) with respect to time. 

2.1. First order equation 

Let us define the microscopic density n{r,t) of a classical many-body system by the re¬ 
lation (ED- The calculation of the sums of type J2s / where / (R^) are any “the one- 
particle” functions, will perform according to the rule 

^/(R^(t))=^y f {r) S{r -Rs{t)) dr = J f (r) n{r,t) dr. (3) 

The first derivative of n (r, t) with respect to time is related to the instantaneous velocities 
of particles. After differentiating the local density n{r,t) with respect to time, we have 

= * e-' (4) 

dt J {2Trf ^ \ dt J 

Using © and Fourier representation of delta-function, we obtain 

■v(R,t)^ n(R,t) dR= -V(v(r,t) n{r,t)). (5) 

This is none other than well known equation of continuity, i.e. a local conservation law of the 
particles number in the system. 

^^ + V(v(r,t)n(r,t))=0. (6) 

Let us pass to derivation of a second order equation. 


2.2. Second order equation 

Let us write a second derivative of the function n(r, t) with respect to time 


d'^n(r,t) 

dt^ 


dk 


^2 k r 


^ (k • R,(t))" - * (k • R,(t)) I . 


( 7 ) 


The first summand in the integrand of this expression can be transformed to the following 
integral 


-E' 


,-ikRa(t) 


k R, 


3 


i(R, t) '?;^(R, t) dR. 


( 8 ) 


S 


3 










To calculate the second summand in the integrand 0, we should evaluate Rs • According 
to the Newton’s second law, we have: 


Rs(t) = — Vr 
m 


= -Vr. 

m 


w (Rs(t) - Rs' {t -t)) + ip (Rs, t) 

. s' 

tT(Rs -R')«(R',<-r(|Rs -R'D) dll' + ip {Rs, t) 


( 9 ) 


where W (Rs — Rs') is an interaction potential between of resting particles located in points Rs 
and Rs', r (|Rs — R'|) is time retardation between points Rg and Rs', (fi (r, t) is an external 
field potential, m is a particle mass, . 

Substituting ® into the second summand in the integrand of 0, we obtain 




I 

m 


e-^^^n{R,t) 


k-VR / lT(R-R')n(R',t-T(|Rs-R'|)) dR' 


( 10 ) 


(k-VR(p{R,t)) 


dR. 


Using ®, 0, and 0 leads to following results 
f dk 


i k 


1 


(27r) 




,-ikRs(t) 


|(k-Rs(t))'| 


1 r d^S (r - R) 
3./ dr'^ 


t(R, t) u^(R, t) dR 


= 3 A [nir,t)v^{r,t)] 

(A is the Laplace operator) and 
f dk 


J 

m 


kr ^g-*kR,(t) |'k.R,(t)^| 

s 

n(r, t) ^Vr ly lU (r — R) n (R, t — t (|Rs — R^|)) dR + ip (r, t) 


( 11 ) 


( 12 ) 


Substituting expressions (HH), (HH) into 0, we obtain the basic equation 

9^n(r,t) 1 * r / N 9/ 

—Q^ = 3 ^ [nir,t)v (r,t)] 


-Vr 


n(r, t) ^Vr Ij W {v — R)n (R, t — r (|r — R|)) dR + (r, t) 


(13) 


The further part of the present paper is devoted to the analysis of this equation and its 
consequences. 
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3. Dynamics of slightly inhomogeneous medium 


Let us consider the simplest limiting case of the equation (1131) — limit a “nearly uniform” 
medium without external field ((/? (r, t) = 0 ): 

/ n(r,t) = no + ni(r,t), no = const, |ni(r, t)| < no, 
t>^(r, t) = const = 

After linearization of the equations m with respect to ni(r,t), we have 

^ = ^Ani(r,t) + —A [ W (r-R)n.i (R, f - t (|r - R|)) dR. (15) 

ot^ 3 rn J 

This linear integro-differential equation is different from the usual hyperbolic wave equation 
by the presence of an integral term of convolution type, so the solution can be found using 
the Fourier integral 


( 271 -) 

Hence, we obtain an equation for hi(k, w) 


, , /■ dkdw _ , , ,, ,, 

/ 7 —^7ii(k,w) exp[i(k-r-a;t)] . 
J ( Stt ) 


(16) 


dk.duj ^ . ,, 

-^ Til (k, (jj) exp [z (k • r — ujt)] 


X ^ 0.2 - ^ 

3 m 


( 2 ^) 
^dRlT(R) 


(17) 


= 0 . 


Thus, the dispersion law oj (k) of oscillations in slightly inhomogeneous medium consisting of 
particles interacting via retarded potentials, satisfies the following equation 


0.2 — — — — fc 2 lTi (k, uj) = 0 , 

3 m 


where 


lTi(k,w) = J dRlT(R) 

is Fourier-transform of complex-valued effective potential 

lFi(r,a;) = lT(r) 


(18) 

(19) 

( 20 ) 


Note that in the absence of interactions delay the dispersion law has the following form 


u) (k) = ±fc 



-W^(k), 


where 


W (k) = IFi (k, w) 




is the real-valued Fourier-transform of the interatomic potential 


( 21 ) 

( 22 ) 


5 









In general case, the roots of the characteristic equation (IT^ with respect to w are complex¬ 
valued functions of the wave vector k. This fact leads to appearance of both the damped 
harmonics and the growing harmonics in (jl6p . Thus, the imaginary parts of the characteristic 
equation roots lead to irreversible behavior of solutions of the equation (HID. 

The most essential qualitative consequence of the equation dT51) is the fact that the func¬ 
tion w(k) is a complex-valued function. This means that retarded interaction between the 
particles leads to irreversible behavior of a classical many-body system. Irreversible determin¬ 
istic behavior of a many-body system takes place without any probabilistic assumptions and 
hypotheses. 


4. The first non-vanishing correction due to the interaction retardation 

In the case of finite range (or rapidly decreasing at large distances) interatomic potentials 
it is enough to consider the first term due to interaction retardation 

u (R, t - T (|r - R|)) ~ n (R, t) - t (|r - R|). (23) 


Substitute this expression into the equation (1151) and obtain 

^ - —A f W{r- R)ni(R,t) dR 

ot‘‘ 3 m J 

+ — A f Wir - R) T (|r - R|) dR = 0. 

m J at 


(24) 


The first two terms on the left side of this equation correspond to the conventional hyper¬ 
bolic wave equation. The third term takes into account the contribution of interatomic forces 
without regard to the retardation. The last term takes into account the retardation and has 
a first derivative with respect to time. The form of this term is similar to a frictional force 
(generally speaking, with alternating friction coefficient). 

Note that equation (I24p is somewhat similar to the well-known hyperbolic heat equation 
(in other words, the heat waves equation) [HQ, which provides a finite speed of heat 
propagation in contrast to the parabolic heat conduction equation. 

Let us represent the solution ni (r, t) of the equation (IMl) in the form of the Fourier integral 

/ dk 

(25) 

y2tlT J 

As a result, we have the following equation with respect to ni(k, t): 


5^ni(k, t) no 


df^ 


_ ■ALk^W2{\^) 


9ni(k, t) 

Wt 


^ + -wik) 

o m 


ni(k,t) = 0, 


(26) 


where ni(k,<), W(k), and W 2 ik) are the Fourier transforms of the functions ni(r,t), IF(r), 
and IT(r)r(|r|), respectively. 
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The roots of the corresponding characteristic equation in the linear approximation with 
respect to r have the form 


= ±^k. 


\ + -w{\.) 

3 m 


no 

.2m 


kW2{k) 


n 2 


^k^W2ik) ±iksj\ + ^fT(k). 

Im Vo m 


(27) 


Depending on the sign of W 2 (k), the corresponding Fourier harmonics are damped oscillations 
(for W 2 (k) < 0), or divergent oscillations (for 1^2 (k) > 0). Thus, even an arbitrarily small 
retardation of the interactions leads to irreversible behavior of many-particle systems. 


5. Conclusion 


The paper contains the following results. 


• A new form of the equations of motion of many-particle system in classical mechanics 
is obtained. Within this approach the desired function is the exact microscopic density 
as a function of position and time. The effect of the interactions between particles 
retardation is taken into account. 

• Derivation and interpretation of obtained equations of motion does not use any statistical 
and probabilistic assumptions. 

• It is shown that the retardation of inter-particle interactions in the many-body system 
leads to irreversibility of the system. 

The basic equation (I13|) is not closed, because it contains two unknown functions: the 
microscopic density n{r,t) and velocity distribution function v(r,t). The second of these 
functions under certain conditions can be expressed via the velocity potential V’(r, f) 

v(r,t) = -V'0(r,t). (28) 


In this case, a closed system of equations of motion for many-body system (in external fields 
absence) has the following form 


dn{r,t) 
dt 


- V [n(r,t) VV'(r,t)] = 0; 


d'^n(r,t) I^ 

dt^ . 3 


Vr 

m 


n{v,t) {\7'ip{r,t)) 

n(r, t) ^Vr |y IF (r — R) n (R, t — r (|r — R|)) c?R 


(29) 


= 0 . 


In general, this system of equations is the subject to study. The continuity equation is invariant 
under the transformation t —>■ —t. The basic equation is not invariant with respect to this 
transformation. Therefore, it can be argued that the presence of interactions delay leads to 
irreversible behavior of many-body systems. 

In this context, there is a lot of interesting open questions. 
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1. Is the Liouville theorem valid for systems with retarded interactions? 

2. Can the equations of continuum mechanics be derived from Can be derived equations 
of continuum mechanics (for example, the Euler and Navier-Stokes equations) from the 
the basic equation (USD? 

3. Are the variational principles valid for the equations with retarded interactions? 

4. How to derive the equilibrium properties of many-particle systems? It is clear that 
as t —>■ oo all the derivatives with respect to time vanish in the basic equation (1131) . Then 
the constant v'^ should be replaced by {k is the Boltzmann constant, T is absolute 
temperature), it follows from elementary kinetic theory. This transformation yields a 
nonlinear integral-differential equation with respect to the microscopic density [3 El. 
But still not quite clear how to calculate the full set of the thermodynamic properties 
of the many-body systems. 
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